Nonsystematic Reporting Biases of the SARS-CoV-2 Variant Mu Could Impact Our Understanding of the Epidemiological Dynamics of Emerging Variants.

Developing a timely and effective response to emerging SARS-CoV-2 variants of concern (VOCs) is of paramount public health importance. Global health surveillance does not rely on genomic data alone to identify concerning variants when they emerge. Instead, methods that utilize genomic data to estimate the epidemiological dynamics of emerging lineages have the potential to serve as an early warning system. However, these methods assume that genomic data are uniformly reported across circulating lineages. In this study, we analyze differences in reporting delays among SARS-CoV-2 VOCs as a plausible explanation for the timing of the global response to the former VOC Mu. Mu likely emerged in South America in mid-2020, where its circulation was largely confined. In this study, we demonstrate that Mu was designated as a VOC ∼1 year after it emerged and find that the reporting of genomic data for Mu differed significantly than that of other VOCs within countries, states, and individual laboratories. Our findings suggest that nonsystematic biases in the reporting of genomic data may have delayed the global response to Mu. Until they are resolved, the surveillance gaps that affected the global response to Mu could impede the rapid and accurate assessment of future emerging variants.

Incorporating variant frequencies data into short-term forecasting for COVID-19 cases and deaths in the USA: a deep learning approach.

BACKGROUND: Since the US reported its first COVID-19 case on January 21, 2020, the science community has been applying various techniques to forecast incident cases and deaths. To date, providing an accurate and robust forecast at a high spatial resolution has proved challenging, even in the short term. METHOD: Here we present a novel multi-stage deep learning model to forecast the number of COVID-19 cases and deaths for each US state at a weekly level for a forecast horizon of 1-4 weeks. The model is heavily data driven, and relies on epidemiological, mobility, survey, climate, demographic, and SARS-CoV-2 variant frequencies data. We implement a rigorous and robust evaluation of our model-specifically we report on weekly performance over a one-year period based on multiple error metrics, and explicitly assess how our model performance varies over space, chronological time, and different outbreak phases. FINDINGS: The proposed model is shown to consistently outperform the CDC ensemble model for all evaluation metrics in multiple spatiotemporal settings, especially for the longer-term (3 and 4 weeks ahead) forecast horizon. Our case study also highlights the potential value of variant frequencies data for use in short-term forecasting to identify forthcoming surges driven by new variants. INTERPRETATION: Based on our findings, the proposed forecasting framework improves upon the available state-of-the-art forecasting tools currently used to support public health decision making with respect to COVID-19 risk. FUNDING: This work was funded the NSF Rapid Response Research (RAPID) grant Award ID 2108526 and the CDC Contract #75D30120C09570.

Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection

The chronic infection hypothesis for novel SARS-CoV-2 variant emergence is increasingly gaining credence following the appearance of Omicron. Here we investigate intrahost evolution and genetic diversity of lineage B.1.517 during a SARS-CoV-2 chronic infection lasting for 471 days (and still ongoing) with consistently recovered infectious virus and high viral genome copies. During the infection, we find an accelerated virus evolutionary rate translating to 35 nucleotide substitutions per year, approximately two-fold higher than the global SARS-CoV-2 evolutionary rate. This intrahost evolution result in the emergence and persistence of at least three genetically distinct genotypes suggesting the establishment of spatially structured viral populations continually reseeding different genotypes into the nasopharynx. Finally, we track the temporal dynamics of genetic diversity to identify advantageous mutations and highlight hallmark changes for chronic infection. Our findings demonstrate that untreated chronic infections accelerate SARS-CoV-2 evolution, providing an opportunity for the emergence of genetically divergent variants. Graphical To understand the intrahost evolution of SARS-CoV-2 from a single patient chronically infected for at least 471 days, Chaguza et al. use whole genome sequencing to estimate the evolutionary rate, the genetic divergence of viral lineages, relative mutation rates, and frequency of mutational variants during the course of the infection.

Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection.

The chronic infection hypothesis for novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant emergence is increasingly gaining credence following the appearance of Omicron. Here, we investigate intrahost evolution and genetic diversity of lineage B.1.517 during a SARS-CoV-2 chronic infection lasting for 471 days (and still ongoing) with consistently recovered infectious virus and high viral genome copies. During the infection, we find an accelerated virus evolutionary rate translating to 35 nucleotide substitutions per year, approximately 2-fold higher than the global SARS-CoV-2 evolutionary rate. This intrahost evolution results in the emergence and persistence of at least three genetically distinct genotypes, suggesting the establishment of spatially structured viral populations continually reseeding different genotypes into the nasopharynx. Finally, we track the temporal dynamics of genetic diversity to identify advantageous mutations and highlight hallmark changes for chronic infection. Our findings demonstrate that untreated chronic infections accelerate SARS-CoV-2 evolution, providing an opportunity for the emergence of genetically divergent variants.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Sequencing SARS-CoV-2 genomes from saliva.

Genomic sequencing is crucial to understanding the epidemiology and evolution of Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Often, genomic studies rely on remnant diagnostic material, typically nasopharyngeal (NP) swabs, as input into whole-genome SARS-CoV-2 next-generation sequencing pipelines. Saliva has proven to be a safe and stable specimen for the detection of SARS-CoV-2 RNA via traditional diagnostic assays; however, saliva is not commonly used for SARS-CoV-2 sequencing. Using the ARTIC Network amplicon-generation approach with sequencing on the Oxford Nanopore MinION, we demonstrate that sequencing SARS-CoV-2 from saliva produces genomes comparable to those from NP swabs, and that RNA extraction is necessary to generate complete genomes from saliva. In this study, we show that saliva is a useful specimen type for genomic studies of SARS-CoV-2.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 variants Alpha and Iota.

SARS-CoV-2 variants shaped the second year of the COVID-19 pandemic and the discourse around effective control measures. Evaluating the threat posed by a new variant is essential for adapting response efforts when community transmission is detected. In this study, we compare the dynamics of two variants, Alpha and Iota, by integrating genomic surveillance data to estimate the effective reproduction number (Rt) of the variants. We use Connecticut, United States, in which Alpha and Iota co-circulated in 2021. We find that the Rt of these variants were up to 50% larger than that of other variants. We then use phylogeography to show that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of Alpha were larger than those resulting from Iota introductions. By monitoring the dynamics of individual variants throughout our study period, we demonstrate the importance of routine surveillance in the response to COVID-19.

COVID-19 Outcomes and Genomic Characterization of SARS-CoV-2 Isolated From Veterans in New England States: Retrospective Analysis.

BACKGROUND: Clinical and virologic characteristics of COVID-19 infections in veterans in New England have not been described. The average US veteran is a male older than the general US population. SARS-CoV-2 infection is known to cause poorer outcomes among men and older adults, making the veteran population an especially vulnerable group for COVID-19. OBJECTIVE: This study aims to evaluate clinical and virologic factors impacting COVID-19 outcomes. METHODS: This retrospective chart review included 476 veterans in six New England states with confirmed SARS-CoV-2 infection between April and September 2020. Whole genome sequencing was performed on SARS-CoV-2 RNA isolated from these veterans, and the correlation of genomic data to clinical outcomes was evaluated. Clinical and demographic variables were collected by manual chart review and were correlated to the end points of peak disease severity (based on oxygenation requirements), hospitalization, and mortality using multivariate regression analyses. RESULTS: Of 476 veterans, 274 had complete and accessible charts. Of the 274 veterans, 92.7% (n=254) were men and 83.2% (n=228) were White, and the mean age was 63 years. In the multivariate regression, significant predictors of hospitalization (C statistic 0.75) were age (odds ratio [OR] 1.05, 95% CI 1.03-1.08) and non-White race (OR 2.39, 95% CI 1.13-5.01). Peak severity (C statistic 0.70) also varied by age (OR 1.07, 95% CI 1.03-1.11) and O2 requirement on admission (OR 45.7, 95% CI 18.79-111). Mortality (C statistic 0.87) was predicted by age (OR 1.06, 95% CI 1.01-1.11), dementia (OR 3.44, 95% CI 1.07-11.1), and O2 requirement on admission (OR 6.74, 95% CI 1.74-26.1). Most (291/299, 97.3%) of our samples were dominated by the spike protein D614G substitution and were from SARS-CoV-2 B.1 lineage or one of 37 different B.1 sublineages, with none representing more than 8.7% (26/299) of the cases. CONCLUSIONS: In a cohort of veterans from the six New England states with a mean age of 63 years and a high comorbidity burden, age was the largest predictor of hospitalization, peak disease severity, and mortality. Non-White veterans were more likely to be hospitalized, and patients who required oxygen on admission were more likely to have severe disease and higher rates of mortality. Multiple SARS-CoV-2 lineages were distributed in patients in New England early in the COVID-19 era, mostly related to viruses from New York State with D614G mutation.

Evaluation of saliva self-collection devices for SARS-CoV-2 diagnostics.

BACKGROUND: There is an urgent need to expand testing for SARS-CoV-2 and other respiratory pathogens as the global community struggles to control the COVID-19 pandemic. Current diagnostic methods can be affected by supply chain bottlenecks and require the assistance of medical professionals, impeding the implementation of large-scale testing. Self-collection of saliva may solve these problems, as it can be completed without specialized training and uses generic materials. METHODS: We observed 30 individuals who self-collected saliva using four different collection devices and analyzed their feedback. Two of these devices, a funnel and bulb pipette, were used to evaluate at-home saliva collection by 60 individuals. SARS-CoV-2-spiked saliva samples were subjected to temperature cycles designed to simulate the conditions the samples might be exposed to during the summer and winter seasons and sensitivity of detection was evaluated. RESULTS: All devices enabled the safe, unsupervised self-collection of saliva. The quantity and quality of the samples received were acceptable for SARS-CoV-2 diagnostic testing, as determined by human RNase P detection. There was no significant difference in SARS-CoV-2 nucleocapsid gene (N1) detection between the freshly spiked samples and those incubated with the summer and winter profiles. CONCLUSION: We demonstrate inexpensive, generic, buffer free collection devices suitable for unsupervised and home saliva self-collection.

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA.

The SARS-CoV-2 Delta variant rose to dominance in mid-2021, likely propelled by an estimated 40%-80% increased transmissibility over Alpha. To investigate if this ostensible difference in transmissibility is uniform across populations, we partner with public health programs from all six states in New England in the United States. We compare logistic growth rates during each variant's respective emergence period, finding that Delta emerged 1.37-2.63 times faster than Alpha (range across states). We compute variant-specific effective reproductive numbers, estimating that Delta is 63%-167% more transmissible than Alpha (range across states). Finally, we estimate that Delta infections generate on average 6.2 (95% CI 3.1-10.9) times more viral RNA copies per milliliter than Alpha infections during their respective emergence. Overall, our evidence suggests that Delta's enhanced transmissibility can be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on underlying population attributes and sequencing data availability.

Longitudinal Immune Profiling of a Severe Acute Respiratory Syndrome Coronavirus 2 Reinfection in a Solid Organ Transplant Recipient.

BACKGROUND: The underlying immunologic deficiencies enabling severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reinfection are currently unknown. We describe deep longitudinal immune profiling of a transplant recipient hospitalized twice for coronavirus disease 2019 (COVID-19). METHODS: A 66-year-old male renal transplant recipient was hospitalized with COVID-19 March 2020 then readmitted to the hospital with COVID-19 233 days after initial diagnosis. Virologic and immunologic investigations were performed on samples from the primary and secondary infections. RESULTS: Whole viral genome sequencing and phylogenetic analysis revealed that viruses causing both infections were caused by distinct genetic lineages without evidence of immune escape mutations. Longitudinal comparison of cellular and humoral responses during primary SARS-CoV-2 infection revealed that this patient responded to the primary infection with low neutralization titer anti-SARS-CoV-2 antibodies that were likely present at the time of reinfection. CONCLUSIONS: The development of neutralizing antibodies and humoral memory responses in this patient failed to confer protection against reinfection, suggesting that they were below a neutralizing titer threshold or that additional factors may be required for efficient prevention of SARS-CoV-2 reinfection. Development of poorly neutralizing antibodies may have been due to profound and relatively specific reduction in naive CD4 T-cell pools. Seropositivity alone may not be a perfect correlate of protection in immunocompromised patients.

Longitudinal immune profiling of a SARS-CoV-2 reinfection in a solid organ transplant recipient.

The underlying immunologic deficiencies enabling SARS-CoV-2 reinfections are currently unknown. Here we describe a renal-transplant recipient who developed recurrent, symptomatic SARS-CoV-2 infection 7 months after primary infection. To elucidate the immunological mechanisms responsible for reinfection, we performed longitudinal profiling of cellular and humoral responses during both primary and recurrent SARS-CoV-2 infection. We found that the patient responded to the primary infection with transient, poor-quality adaptive immune responses that was further compromised by intervening treatment for acute rejection of the renal allograft prior to reinfection. Importantly, we identified the development of neutralizing antibodies and humoral memory responses prior to SARS-CoV-2 reinfection. However, these neutralizing antibodies failed to confer protection against reinfection, suggesting that additional factors are required for efficient prevention of SARS-CoV-2 reinfection. Further, we found no evidence supporting viral evasion of primary adaptive immune responses, suggesting that susceptibility to reinfection may be determined by host factors rather than pathogen adaptation.

Diverse functional autoantibodies in patients with COVID-19.

COVID-19 manifests with a wide spectrum of clinical phenotypes that are characterized by exaggerated and misdirected host immune responses1-6. Although pathological innate immune activation is well-documented in severe disease1, the effect of autoantibodies on disease progression is less well-defined. Here we use a high-throughput autoantibody discovery technique known as rapid extracellular antigen profiling7 to screen a cohort of 194 individuals infected with SARS-CoV-2, comprising 172 patients with COVID-19 and 22 healthcare workers with mild disease or asymptomatic infection, for autoantibodies against 2,770 extracellular and secreted proteins (members of the exoproteome). We found that patients with COVID-19 exhibit marked increases in autoantibody reactivities as compared to uninfected individuals, and show a high prevalence of autoantibodies against immunomodulatory proteins (including cytokines, chemokines, complement components and cell-surface proteins). We established that these autoantibodies perturb immune function and impair virological control by inhibiting immunoreceptor signalling and by altering peripheral immune cell composition, and found that mouse surrogates of these autoantibodies increase disease severity in a mouse model of SARS-CoV-2 infection. Our analysis of autoantibodies against tissue-associated antigens revealed associations with specific clinical characteristics. Our findings suggest a pathological role for exoproteome-directed autoantibodies in COVID-19, with diverse effects on immune functionality and associations with clinical outcomes.

Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2.

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69-70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1.

Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Detection of SARS-CoV-2 RNA by multiplex RT-qPCR.

The current quantitative reverse transcription PCR (RT-qPCR) assay recommended for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) testing in the United States requires analysis of 3 genomic targets per sample: 2 viral and 1 host. To simplify testing and reduce the volume of required reagents, we devised a multiplex RT-qPCR assay to detect SARS-CoV-2 in a single reaction. We used existing N1, N2, and RP primer and probe sets by the Centers for Disease Control and Prevention, but substituted fluorophores to allow multiplexing of the assay. The cycle threshold (Ct) values of our multiplex RT-qPCR were comparable to those obtained by the single assay adapted for research purposes. Low copy numbers (≥500 copies/reaction) of SARS-CoV-2 RNA were consistently detected by the multiplex RT-qPCR. Our novel multiplex RT-qPCR improves upon current single diagnostics by saving reagents, costs, time, and labor.

Real-time public health communication of local SARS-CoV-2 genomic epidemiology.

Genomic epidemiology can provide a unique, real-time understanding of SARS-CoV-2 transmission patterns. Yet the potential for genomic analyses to guide local policy and community-based behavioral decisions is limited because they are often oriented towards specially trained scientists and conducted on a national or global scale. Here, we propose a new paradigm: Phylogenetic analyses performed on a local level (municipal, county, or state), with results communicated in a clear, timely, and actionable manner to strengthen public health responses. We believe that presenting results rapidly, and tailored to a non-expert audience, can serve as a template for effective public health response to COVID-19 and other emerging viral diseases.